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S ummar y 

This paper presents comparative analyses of atmospheric radon (Rn) distributions simulated 
using different emission scenarios and the observations. Results indicate that the model generally 
reproduces observed distributions of Rn but there are some biases in the model related to 
differences in large-scale and convective transport. Simulations presented here use an off-line 
three-dimensional chemical transport model driven by assimilated winds and two scenarios of Rn 
fluxes (atom cm'V 1 ) from ice- free land surfaces: (A) globally uniform flux of 1.0, and (B) 
uniform flux of 1.0 between 60°S and 30°N followed by a sharp linear decrease to 0.2 at 70°N. 
We considered an additional scenario (C) where Rn emissions for case A were uniformly 
reduced by 28%. Results show that case A overpredicts observed Rn distributions in both 
hemispheres. Simulated northern hemispheric (NH) Rn distributions from cases B and C 
compare better with the observations, but are not discernible from each other. In the southern 
hemisphere, surface Rn distributions from case C compare better with the observations. We 
performed a synoptic scale source-receptor analysis for surface Rn to locate regions with ratios 
B/A and B/C less than 0.5. Considering an uncertainty in regional Rn emissions of a factor of 
two, our analysis indicates that additional measurements of surface Rn particularly during April- 
October and north of 50°N over the Pacific as well as Atlantic regions would make it possible to 
determine if the proposed latitude gradient in Rn emissions is superior to a uniform flux scenario. 
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ABSTRACT 


We present comparative analyses of atmospheric radon (Rn) distributions simulated using 
different emission scenarios and the observations. Results indicate that the model generally 
reproduces observed distributions of Rn but there are some biases in the model related to 
differences in large-scale and convective transport. Simulations presented here use an off-line 
three-dimensional chemical transport model driven by assimilated winds and two scenarios of Rn 
fluxes (atom cm V 1 ) from ice-free land surfaces: (A) globally uniform flux of 1.0, and (B) 
uniform flux of 1.0 between 60°S and 30°N followed by a sharp linear decrease to 0.2 at 70°N. 
We considered an additional scenario (C) where Rn emissions for case A were uniformly 
reduced by 28%. Results show that case A overpredicts observed Rn distributions in both 
hemispheres. Simulated northern hemispheric (NH) Rn distributions from cases B and C 
compare better with the observations, but are not discernible from each other. In the southern 
hemisphere, surface Rn distributions from case C compare better with the observations. We 
performed a synoptic scale source-receptor analysis for surface Rn to locate regions with ratios 
B/A and B/C less than 0.5. Considering an uncertainty in regional Rn emissions of a factor of 
two, our analysis indicates that additional measurements of surface Rn particularly during April- 
October and north of 50°N over the Pacific as well as Atlantic regions would make it possible to 
determine if the proposed latitude gradient in Rn emissions is superior to a uniform flux scenario. 
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INTRODUCTION 


Evaluation of the representation of large- and sub-scale transport processes in an 
atmospheric model is a prerequisite for its application to study chemistry-climate interactions. 
Long lived tracers such as CFCs and SF6 are generally used to evaluate large-scale atmospheric 
transport features such as interhemispheric exchange (Prather et al., 1987; Denning et al., 1999; 
Gupta et al., 2001). Radon ( Rn, or simply Rn), characterized by its short radioactive decay 
time constant (ti/2=3.825 days) and primarily land-based emissions, has been extensively used as 
a tracer for study of sub-scale convective and synoptic transport processes in the troposphere 
(Liu et al., 1984; Dentener et al., 1999; Chevillard et al., 2002). Radon has been used for model 
intercomparison studies (Jacob et al., 1997; Rasch et al., 2000), assessment of convective 
parameterizations (Mahowald et al., 1995; Allen et al., 1996), and evaluation of model transport 
of continental air to remote oceanic locations (Heimann et al., 1990; Balkanski and Jacob, 1990). 

Emissions of Rn from radioactive decay of radium in soils are known to be highly 
variable. Factors controlling Rn emissions include the 226 Ra concentration and the physical 
characteristics of soil affecting its escape into the atmosphere (Turekian et al., 1977; Nazaroff, 
1992). Jacob and Prather (1990) and Genthon and Armengaud (1995) have derived algorithms to 
account for some of the variations in Rn emissions. Observations used to develop these 
algorithms are limited and do not uniformly represent the entire globe, therefore a simple 
emission rate of 1.0 atom cm'V 1 over the ice-free land surfaces between 60°S and 60°N has been 
widely used (Jacob et al., 1997; Rasch et al., 2000]. This emission rate is considered to be 
accurate to within 25% globally and within a factor of two regionally (see Jacob et al., 1997 and 
references therein). Also, this Rn flux falls between two other estimates of 0.72 atom cm'V 1 and 
1.2 atom cm'V 1 , suggested by Lambert et al. (1982) and Turekian et al. (1977). Various models 
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have simulated atmospheric distributions of Rn using this emission rate and assimilated 
meteorology and compared them against observed surface concentrations and vertical profiles 
(e.g., Dentener et al., 1999 and Taguchi et al., 2002). Oceanic contributions to the Rn flux are 
considered negligibly small (0.005 atom cm'V 1 ). 

Recently, Conen and Robertson (2002) reported that land-based Rn fluxes decrease 
linearly from 1.0 atom cm'V 1 at 30°N to 0.2 atom cm'V 1 at 70°N. Their conclusion was based 
on indirect estimates of Rn emissions. The authors speculated that the northward decrease in Rn 
emissions could be due to more organic soils with decreasing radium content and shallow water 
tables. They emphasized the potential importance of this finding in validation of atmospheric 
transport models. This revised estimate modifies Rn emission from 59% of the northern 
hemispheric (NH) land-cover between 31°N and 71°N. 

Here, we investigate the impact of different emission estimates of Rn on simulated 
atmospheric distributions. We compare simulated Rn distributions against corresponding 
observations. This comparative analysis will also evaluate treatment of physical processes and 
meteorological inputs to the model that are responsible for long-range as well as vertical 
transport of Rn in the atmosphere. 

MODEL SIMULATIONS 

All the simulations presented here were performed using the Goddard Space Flight 
Center (GSFC) off-line three-dimensional parameterized chemistry and transport model (PCTM) 
(Nielsen and Douglass, 2001). Meteorological fields (surface pressure, temperature and winds), 
and convective and diffusive fluxes input to PCTM were taken from a reanalysis run of the 
Goddard Earth Observation System, version 4 (GEOS-4) data assimilation system given at a 
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resolution of 2.5°x2° with 55 hybrid sigma-pressure levels between the surface and 80 km. 
GEOS-4 is an updated version of the Finite-Volume Data Assimilation System (FVDAS) used 
by Douglass et al. (2003). The time dependent tracer mass continuity equation was solved at 15- 
minute time intervals using linearly interpolated wind fields drawn from six-hour archives of 
GEOS-4 meteorology. Advection was calculated using a flux-form semi-Lagrangian scheme 
(Lin and Rood, 1996) on the GEOS-4 grid. A pressure fixer, applied following Prather et al. 
(1987), adjusts the horizontal winds slightly to remove the mass increase that results from 
inconsistency between the forecast and observed surface pressures used in the assimilation 
process. A semi implicit numerical scheme was applied for convective tracer transport based on 
the three-dimensional cloud mass flux produced by the FVDAS. 

Atmospheric Rn simulations were performed using analyzed meteorology for 1993 and 
1994. For both years, twoemi ssion scenarios were used: (A) 1.0 atom cm'V 1 (Jacob et al., 1997) 
and O latitude dependent emissions suggested by Conen and Robertson (2002). Fractional land- 
based emissions were uniformly applied in coastal and island-containing model grids using a 
high resolution (l/6°xl/6°) surface topography. No alterations in land-based Rn fluxes in these 
model grids were introduced even though there is an evidence of lower Rn emanation rates at 
island of Hawaii due to water saturation (Whitdestone et al., 1996). Also, surface Rn emissions 
were zero for model grids with surface temperature below freezing point. We lack knowledge of 
the temporal variations in land-based Rn sources, therefore, constant emissions were applied at 
each time step. For both cases, the model was allowed to spin-up for three months before 
comparing the results with observations. We also considered another Rn emission rate of 0.72 
atom cm'V 1 , as proposed by Lambert et al. (1982), by proportionally scaling down Rn 
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concentrations from case A. This model scenario is termed as C. For cases A, B and C, annual 
global Rn emissions were 62.1, 55.3 and 44.7 moles respectively. 

To investigate the impact of latitudinal gradient in Rn emissions, we concentrate on its 
NH distributions. Validation was performed by comparison with observations from the same 
period as the analyzed wind [1993-94]. Measured Rn vertical profiles over Moffett Field during 
1994 (Kritz et al., 1998) and from the North Atlantic Regional Experiment (NARE) Intensive- 
1993 (Zaucker et al., 1996) were compared with model profiles for all three emission scenarios. 
Surface measurements of Rn over Bermuda (64.45°W, 32.2°N, 28m asl; BER) for 1993 reported 
by Environmental Measurement Laboratory (EML) (flutter et al., 1995) were also used. EML 
reports data for BER at a 60-minute time interval. Even though Rn data for other northern 
hemispheric stations e.g. island site (Mauna Loa) and continental sites (Freiburg and 
Schauinsland, Germany) are also available, we considered Bermuda as a representative site for 
model-data comparison because of its remote location, low surface elevation and negligible local 
emissions (surface land area= 53 km 2 ). For the purpose of model evaluation to examine the 
influence of uniform reduction in surface emissions, we also include comparison between 
simulated and observed surface Rn for 1993 at two southern hemispheric (SH) stations: 
Amsterdam Island [77.34°E, 37.5°S, 50m asl; AMS] and Cape Grim [144.41°E, 40.41°S, 94m 
asl; CGA] following Dentener et al. (1999). CGA is located on the coastal site of Australia and 
AMS is remotely located in the Indian Ocean with land surface land area of 61 km 2 . 

RESULTS AND DISCUSSION 
A. Comparison for Rn vertical profiles 

a. Vertical profiles during the NARE- 199 3 campaign 
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Simulated and observed Rn for the NARE Intensive period 16-31 August, 1993 are 
compared in Fig. 1. Sixty-six Rn samples were collected between the boundary layer and 5.5 km 
over the North Atlantic Ocean and the continent in the vicinity of Nova Scotia and New 
Brunswick, Canada (Zaucker et al., 1996). To facilitate the comparison, observed Rn 
concentrations and corresponding modeled values were divided into three subdivisions: 0-2 km 
for the marine boundary layer (MBL); 0-2 km for continental boundary layer (BL); and free 
troposphere (FT). The total numbers of samples in the MBL, BL and FT subdivisions are 28, 13 
and 25 respectively. Average concentrations and ±1 ct standard deviation for each division were 
computed for analysis. 

Case A predicts larger than observed Rn concentrations for all three regions. Modeled Rn 
concentrations from cases B and C closely correspond to observed values, particularly for MBL 
and FT. For case B, the agreement is within ±5% and for case C, it is within -3 and 16%. For the 
BL, cases B and C overpredict observed average Rn by 27% and 38% respectively. This stated 
discrepancy could be due to uncertainty in Rn emissions and insufficient vertical transport of 
continental boundary layer airmass to the free troposphere. The latitudinal band covered during 
the NARE- 1993 Intensive is 42°N-46°N where the Rn flux was reduced by 24-32% in case B. 
The mean value of this reduction corresponds to the decrease in global Rn surface emissions for 
case C as compared to that for case A. This comparison indicates that although the ratio (B/C) 
for integrated NH Rn source strength is 1.25, the simulated Rn distributions from both cases for 
the NARE geographic area are not clearly distinct from each other. 

b. Vertical profiles over Moffett Field during 1994 
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A comparison between modeled and observed Rn concentrations over Moffett Field 
(37.4°N, 122.0^^0 is shown in Fig. 2. Observations are adopted from Kritz et al. (1998), who 
report Rn profiles from the surface to 11.5 km obtained from 127 samples on 11 flights during 
June 3-August 16, 1994. We averaged Rn concentrations corresponding to three segments of 
vertical profiles: (1) 0-2 km; (2) 2-6 km; and (3) 6-11.5 km. All three emission scenarios 
overestimate observed average Rn in the lowest 2 km of the atmosphere, however, this 
comparison is the worst for case A. In the model, Moffett Field is located in the coastal grid 
(with 64% of land cover) which, as stated earlier, may have lower Rn emanation rate than what 
is considered here. At the latitude of Moffett Field, reductions in Rn emission for cases B and C 
are 15% and 28% respectively relative to case A. For the layer between 2-6 km, Rn distributions 
from all cases closely follow each other within ±la deviation around the mean. Besides 
atmospheric mixing, a combination of two factors: (1) higher magnitude of NH emissions, and 
(2) a relatively smaller reduction in local surface flux at Moffett field for case B resulted in its 
distributions thata re not different from those of case C. 

Because of simulated excessive Rn near the surface, irrespective of emission scenarios, 
model underpredicts observed Rn at high altitudes. This discrepancy may be due to insufficient 
long-range transport and possible existence of high Rn emission from East Asian region 
(Dentener et al., 1999) which is not considered here. Doubling of convective fluxes in model 
simulation improved the results by up to 15% at high altitudes. Model simulations for case B 
with regionally tagged emissions for Indo-China and North American regions show that, above 6 
km, more than 50% of the Rn over Moffett Field originates in the Indo-China region (Fig. 3). 
This result is qualitatively consistent with the conclusions of Kritz et al. (1990) and Stockwell et 
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al. (1998). Togetherthe se regional sources account for more than 80% of total atmospheric Rn 
simulated over the Moffett Field. 

B. Comparison for surface Radon 

a. Northern Hemisphere 

We compare modeled and observed Rn surface concentrations for NH station BER. 
Bermuda is downwind of North America and routinely experiences strong influence of its Rn 
sources due to westerly flow. As should be expected due to use of assimilateckneteorological 
fields, the model often captures observed synoptic scale variability in surface Rn as shown for 
June in Fig. 4. In this figure, we used model results for case B, however, it should be emphasized 
that synoptic variability in simulated Rn is independent of choice of emission scenario. However, 
due to differences in emission intensities, the simulated extreme Rn concentrations during these 
synoptic events differ significantly. The model failed to simulate a prolonged event of high Rn 
concentration observed around Julian days 169 and 171. It also predicts excessve Rn around 
Julian days 160-163, but generally corresponds with a data within a factor of two. Radon 
simulations performed using assimilated winds from an earlier version of GEOS exhibited 
similar nature of synoptic scale variability at this location (Alen et al. , 1996). 

For comparison of modeled and observation data for BER over a year, surface Rn was 
averaged each month. The model produces the general features of the observed monthly 
variation in surface Rn at BER (Fig. 5). This figure also displays ±lc deviation about the mean 
of the observations. Standard error for the observation is calculated to be very small. Lower Rn 
concentrations during late summer and fall are due to the formation of the so-called ‘Bermuda 
High’ which prevents outside air from entering into the high pressure region. 
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During the first half of the year, case A consistently oveipredicts observed monthly 
averaged surface Rn at BER. However, its monthly mean values liewithin the + lo deviation 
around the observed mean values, except for April. Simulated surface Rn concentrations from 
cases B and C closely follow each other and are in close agreement with the observed values 
(with correlation coefficients of 0.60 and 0.55 respectively) despite a reduction of only 6% in 
surface flux at 33°N for case B as compared to 28% for case C. It should be noted these 
reductions in Rn emissions are within the range of uncertainties in regional and global sources. 
However, in the absolute sense, this comparison reinforces previous conclusions that despite a 
difference in NH source strength in Rn emissions for cases B and C, their simulated Rn 
distributions do not significantly differ from each other due to mixing of airmass from other 
latitudes. 

b. Southern hemisphere 

In this section, we present comparison between simulated and observed surface Rn for 
the two southern hemispheric stations, AMS and CGA. Fig. 6 shows a comparison of synoptic 
changes in observed and simulated surface Rn from case B at (a) CGA for Marchand (b) AMS 
for July. Similar to that for BER station, at both sites the model captures most events of observed 
high and low surface Rn activities. In case of AMS, the model successfully simulates the 
baseline Rn surface concentrations and events of radonic storms. In the absence of local 
emissions due to very small surface land area, the good comparison of these characteristic 
features of surface Rn at AMS establishes a confidence in model’s ability to accurately simulate 
the transport of continental airmass to the remote sites. 

Figure 7 shows a comparison between simulated and observed monthly averaged surface 
Rn at (a) CGA and (b) AMS. This figure also displays ±la deviation about the mean for 
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observation. Standard error for the observation is calculated to be very small. At CGA, the model 
reproduces monthly variation in observed surface Rn within ±1 standard deviation from the 
mean except for the October month. For this month both cases A and B overpredict the surface 
Rn. In general, simulated monthly mean Rn from case C closely follows corresponding observed 
surface values at CGA. For case C, coefficient of correlation between monthly averaged 
observed and simulated surface Rn is 0.68. There seems to be a tendency in the model to 
underpredict the observed mean Rn during May and June when corresponding values from all 
three cases are significantly lower. 

At AMS, all cases reproduce the general features of monthly variation in observed 
surface Rn but overpredict its magnitude indicating insufficient ventilation of boundary layer 
airmass to free troposphere. This overprediction is the least for case C but it is within +1 standard 
deviation from the observed mean except for January and December. For this case, the 
coefficient of correlation between monthly averaged observed and simulated surface Rn is 0.72. 
For most of the months, the simulated monthly mean surface Rn from cases A and B lies outside 
one standard deviation from the corresponding observed mean value (Fig. 7b). Combined with 
the inference from the earlier discussion for CGA station, this comparison shows that emission 
scenario from case C is better suited for simulation of SH Rn distribution. 

We performed another model simulation identical to case A except that the effect of Rn 
emissions due to partial land cover in coastal and island grid boxes was neglected. For the AMS 
site, there is no significant effect of neglecting the local emissions due to its small land area as 
shown by plus symbols in Fig. 7b. On the other hand, improper distribution of Rn emissions in 
the model grid containing CGA leads to overestimation of its simulated surface concentrations 
for the months of September-December. Similarly, for two other SH sites, Kerguelen (70.15°E, 
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49.21°S, 30m asl) and Crozet (51.52°E, 46.26°S, 150m asl), results not shown here, model 
calculations indicate that local emissions can significantly affect the baseline surface Rn because 
of their land surface areas (7215 km and 325 km respectively). These results indicate that Rn 
fluxes need to be prescribed carefully in coastal and island containing model grids in order to 
compare to data meaningfully. 

Based on comparisons with the observations presented here, this study clearly shows that the 
model, in general, is capable to reproduce the observed distributions of atmospheric Rn. 
Additionally, comparisons show that the use of either the newly proposed NH latitude-dependent 
Rn flux scenario (case B; Conen and Robertson, 2002) or the uniform global ice-free land-based 
Rn surface flux of 0.72 atom cm V 1 (case C; Lambert et ai., 1982) in model simulations leads to 
NH atmospheric distributions of Rn that are in good agreement with the observations. Because of 
the limited observation database, Rn distributions from these scenarios are not clearly discernible 
from each other. Case A (with continental Rn surface flux 1.0 atom crriV 1 ) overpredicts 
observed Rn concentrations in both hemispheres. At SH sites, Rn distributions from case C 
compare better with the corresponding observations. 

C. Comparison between NH Rn distributions simulated using different emission scenarios 

In this section, we concentrate on how the atmospheric Rn distribution for emission scenario 
B differs from those of A and C in the NH on various time scales. 
a. Annual averaged distributions 

Fig. 8 shows the NH latitude-longitude distribution of the ratio of annually averaged surface 
Rn mixing ratios between cases B and A. In this figure, we also show locations of Rn 
measurement sites that are influenced by the differences in both emission distributions. Because 
of Rn’s short lifetime and the prescribed differences in its NH emissions, values of B/A lower 
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than 1.0 are mostly confined to the NH and decrease towards higher latitudes (shown bjsolid 
lines). South of 45°N, the ratio B/A, in general, is more than 0.75 which is within the range of 
global uncertainty in Rn surface emissions. At high latitudes over continental areas this ratio is 
lower than 0.5. There are two factors that contribute to this pattern of B/A: (1) northward decline 
in Rn emissions, and (2) about 12% less total NH emissions for case B. After uniformly scaling 
down NH land-based emissions for case A by 12%, values of B/A are as low as 0.7 at high 
latitudes (shown by dash lines) indicating that it is not the emission strength but is the latitudinal 
gradient that has strong impact on the distribution of B/A. Calculated lower values of B/A persist 
throughout the troposphere for most of the NH as shown in Fig. 9. Note that the lowest value of 
B/A is centered between 55°N and 60°N. At the mid-latitudes, there is a tendency of an increase 
in B/A values with altitude indicating decline in impact of the latitudinal gradient in Rn 
emissions. At 50°N and 4 km altitude, annually and zonally averaged Rn distributions from 
simulations B and C would not be discernible as both ratios B/A and C/A are roughly equal to 
0.72. 

b. Source-receptor analysis on synoptic timescales 

To further identify regions where synoptic scale variations in surface Rn from case B can be 
consistently and distinedy separated from those of cases A and C, we concentrate on 184 surface 
receptors. These receptors were placed between latitude band of 25°N-73°N and longitude bands 
of 10°W-75°W (Adantic Region), 132.5°E-180°E and 125°W-180°W (Pacific Region) as shown 
in Fig. 10. The size of each receptor is 7.5°x6.0°. We sampled the simulated Rn distributions at 
these receptors from simulations A and B over the entire year at a time resolution of 15 minutes. 
We scaled down Rn distributions from case A at each receptor by 0.72 to generate Rn 
distributions for case C. Considering that regional Rn emissions have an uncertainty of a factor 
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of two, we selected only those synoptic events at each receptor where ratios B/A and B/C of 
surface Rn is less than 0.5. We performed three simulations for case B in which Rn emissions 
were tagged for latitudinal bands: (Bl) 29°N-43°N, (B2) 43°N-57°N and (B3) 57°N-71°N. The 
purpose of these additional simulations is to quantify relative sizes of contributions from these 
latitude bands during the selected synoptic events at these receptors. Here, we concentrate on 
simulations A and B, and B3. Our analysis indicates that during March-September, on an 
average 44% times each receptor placed between 50°N and 60°N in Pacific region experienced 
B/A ratios that were less than 0.5. Most of these selected synoptic events occurred when 
contribution from region B3 to surface Rn from case B was morethan 50% . 

Figure 1 1 shows the percent frequency distribution of selected events each month for (a) 
Pacific region and (b) Atlantic region when the ratio B/A is less than 0.5 and contribution from 
B3 region is more than 50%. Figure 12 displays the similar plots for percent frequency 
distribution for ratio B/C. There are several points that can be noted from these figures. First of 
all, for receptors numbered below 43 (i.e. below 43°N located in the Pacific region), the 
occurrence of the events for ratio B/A constrained by both conditions is very limited over the 
entire year (Fig. 1 la). But, as one moves northwards up to 55°N (i.e. for receptor numbers 43-70) 
this frequency increases dramatically during April-October months. For all these receptors, 
transport of airmass from higher latitudes played a major role in the displayed frequency pattern 
in B/A ratio because difference in local emissions south of 55°N will not solely force this ratio to 
be lessthan 0.5 . There is negligible occurrence of these events all around the year for receptors 
numbered below 149 located in the Atlantic region below 50°N (Fig lib). Above 55°N, both 
conditions are met more frequently for most of the year in Pacific as well as Atlantic regions. As 
should be expected due to lower concentrations of surface Rn for case C relative to those for case 
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A, ratio B/C meets both conditions less frequently for receptors placed above 50°N for both 
regions (Fig 12). Note that most of the receptors between 50°N and 60°N are placed over the 
oceanic surface and are free from uncertainties in in-situ emissions. Based on this analysis and 
NH distribution of Rn emissions below 60°N, additional surface measurements of atmospheric 
Rn between 50°N and 60°N (i.e. for receptors 57-84 and 149-166) over the Pacific and Atlantic 
regions, mainly over the oceanic sites, during April-October would make it possible to confirm 
or refute the latitudinal gradient in Rn source distribution (case B). We further stress that these 
measurements of atmospheric Rn must conform to the same reference to avoid any systematic 
biases as discussed by Colie et al. (1995). 

CONCLUSIONS 

We present results from a study in which simulated atmospheric distributions using three 
scenarios of Rn surface fluxes were compared with each other and with the observations. 
Comparisons indicate that, in general, the model reproduces observed surface and vertical 
profiles of atmospheric Rn but there are some biases in the model related to differences in large- 
scale and convective transport. We show that a globally uniform land-based Rn surface flux of 
1.0 atom cm'V 1 overpredicts atmospheric Rn in both hemispheres. Globally uniform surface flux 
of 0.72 atom cm'V 1 or a scenario in which emission decreases with latitude north of 30°N 
produces atmospheric distributions of Rn that compare well with the observed northern 
hemispheric distributions but are not discernible from each other. However, uniform surface flux 
of 0.72 atom cm'V 1 is a better candidate of Rn emissions in the southern hemisphere. This study 
shows that calibrated measurements of surface Rn between 50°N and 60°N over the Pacific and 
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Atlantic regions during April-October would be able to validate the suggested latitudinal gradient 
in surface Rn emissions. 
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Figure Captions 


Fig. 1. Comparison between simulated and observed Rn profiles during NARE-1993. Error 
bars represent ±lc deviation from the mean. 

Fig. 2. Comparison between simulated and observed Rn vertical profiles over Moffett Field for 
the period June 3- August 16, 1994. Error bars represent ±la deviation from the mean. 

Fig. 3. Percent contribution of regional sources of Rn from North American and Indo-China to 
its vertical profile over Moffett Field during June 3-August 16, 1994. Note that 
combined contribution from both regional sources is more than 80% of the total Rn at 
all vertical levels. 


Fig. 4. Comparison between synoptic scale variation in simulated and observed surface Rn at 
Bermuda for June, 1993. 


Fig. 5. Comparison between simulated and observed monthly averaged surface Rn at Bermuda 
during 1993. Observation data were missing for September and October. Bars represent 
±lc deviation from the observed mean. 


Fig. 6. Comparison between synoptic scale variations in simulated and observed surface Rn at 
(a) Cape Grim (CGA) for March, 1993 and (b) Amsterdam (AMS) for July, 1993. 

Fig. 7. Comparison between simulated and observed monthly averaged surface Rn at (a) Cape 
Grim (CGA) and (b) Amsterdam (AMS) during 1993. Bars represent ±lo deviation 
from the observed mean. Plus symbols correspond to surface Rn distribution from case 
A when fractional land/ocean covers were neglected. 

Fig. 8. Horizontal cross-section of the ratio B/A (solid lines) for annually averaged surface Rn. 
The contour levels decrease northwards with an interval of 0.1. Corresponding ratios 
involving case A in which its NH Rn surface fluxes were uniformly reduced by 12% 
are shown by dash lines. Also shown are locations of some of the NH Rn measurement 
sites affected by differences between emissions for cases A and B. 

Fig. 9. Zonally averaged distribution of ratio B/A for simulated annually averaged Rn. 
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Fig. 10. Geographical locations of 184 regional surface receptors placed within Pacific and 
Atlantic regions between 25°N and 73°N. For each region, receptors are identified by 
their serial number which increases from left to right and from bottom to top. 


Fig. 11. Monthly variations in percent frequency distribution of surface Rn events at receptors 
placed over (a) Pacific and (b) Atlantic regions. Contour levels on these plots are 5, 15, 
25, 50, 75 and 95. These events are characterized by ratio B/A less than 0.5 and more 
than 50% contribution from region B3 (57°N-71°N). Vertical dash lines encloses 
latitude band of 49°N-61°N. Latitude correspondence of surface receptors is shown in 
Fig 10. 


Fig. 12. Same as Fig. 11 but for the percent frequency distribution of occurrence of ratio C/A 
less than 0.5. 
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